arXiv:l507.03414v2 [hep-ex] 25 Jul 2015 


EUROPEAN ORGANIZATION EOR NUCLEAR RESEARCH (CERN) 



CERN-PH-EP-2015-153 
LHCb-PAPER-2015-029 
July 13, 2015 


Observation of J/V’p resonances 
consistent with pentaquark states in 
Al ^ J/^K- 'p decays 

The LHCb collaboratioiJ3 

Abstract 

Observations of exotic structures in the channel, which we refer to as 

charmonium-pentaquark states, in —)• J/'ijjK~p decays are presented. The data 
sample corresponds to an integrated luminosity of 3 fb~^ acquired with the LHCb 
detector from 7 and 8 TeV pp collisions. An amplitude analysis of the three-body 
final-state reproduces the two-body mass and angular distributions. To obtain a 
satisfactory fit of the structures seen in the J/'ij^p mass spectrum, it is necessary 
to include two Breit-Wigner amplitudes that each describe a resonant state. The 
significance of each of these resonances is more than 9 standard deviations. One has 
a mass of 4380 ± 8 ± 29 MeV and a width of 205 ± 18 ± 86 MeV, while the second is 
narrower, with a mass of 4449.8 ± 1.7 ± 2.5 MeV and a width of 39 ± 5 ± 19 MeV. 

The preferred assignments are of opposite parity, with one state having spin 3/2 
and the other 5/2. 
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Introduction and summary 


The prospect of hadrons with more than the minimal qnark content {qq or qqq) was 
proposed by Gell-Mann in 1964 Q and Zweig [^, followed by a quantitative model for two 
quarks plus two antiquarks developed by Jaffe in 1976 The idea was expanded upon 
to include baryons composed of four quarks plus one antiquark; the name pentaquark was 
coined by Lipkin [^. Past claimed observations of pentaquark states have been shown to 
be spurious ^|, although there is at least one viable tetraquark candidate, the Z(4430)’'' 
observed in —)■ 7r~^ decays |7]-[^, implying that the existence of pentaquark baryon 

states would not be surprising. States that decay into charmonium may have particularly 
distinctive signatures (lo) . 

Large yields of —)■ J/'iljK~p decays are available at LHCb and have been used for 
the precise measurement of the lifetime 0 (In this Letter mention of a particular 
mode implies use of its charge conjugate as well.) This decay can proceed by the diagram 
shown in Fig. [^a), and is expected to be dominated by A* —)■ K~p resonances, as are 
evident in our data shown in Fig. [^a). It could also have exotic contributions, as indicated 
by the diagram in Fig. Ilb), that could result in resonant structures in the J/0p mass 
spectrum shown in Fig. ^b). 



Figure 1: Feynman diagrams for (a) A^ —>■ J/tpA* and (b) A^ —)• K decay. 




Figure 2: Invariant mass of (a) K p and (b) J/ipp combinations from A^ —)• J/piK p decays. 
The solid (red) curve is the expectation from phase space. The background has been subtracted. 
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In practice resonances decaying strongly into J/ijjp must have a minimal quark content 
of ccuud, and thus are charmonium-pentaquarks; we label such states P^, irrespective of 
the internal binding mechanism. In order to ascertain if the structures seen in Fig. |^b) 
are resonant in nature and not due to reflections generated by the A* states, it is necessary 
to perform a full amplitude analysis, allowing for interference effects between both decay 
sequences. 

The ht uses hve decay angles and the K~p invariant mass rriKp as independent variables. 
First we tried to ht the data with an amplitude model that contains 14 A* states listed by 
the Particle Data Group 12 . As this did not give a satisfactory description of the data, 
we added one state, and when that was not sufficient we included a second state. The 
two states are found to have masses of 4380 ± 8 ± 29 MeV and 4449.8 ± 1.7 ± 2.5 MeV, 
with corresponding widths of 205 ± 18 ± 86 MeV and 39 ± 5 ± 19 MeV. (Natural units are 
used throughout this Letter. Whenever two uncertainties are quoted the hrst is statistical 
and the second systematic.) The fractions of the total sample due to the lower mass and 
higher mass states are (8.4 ± 0.7 ± 4.2)% and (4.1 ± 0.5 ± 1.1)%, respectively. The best ht 
solution has spin-parity values of (3/2“, 5/2+). Acceptable solutions are also found 
for additional cases with opposite parity, either (3/2+, 5/2“) or (5/2+, 3/2“). The best 
ht projections are shown in Fig. Both rriKp and the peaking structure in mj/^p are 
reproduced by the ht. The signihcances of the lower mass and higher mass states are 9 
and 12 standard deviations, respectively. 


2200 - 
2000 |- 
1800 ^ 
1600 
1400P 
1200 |- 
1000 |- 
800 |- 
600 F 

400 

200R 


(a) LHCb 


. V “ 



-■— data 
-•— total fit 

- background 

^ P,( 4450 ) 
^ Pc( 4380 ) 
■+■■ A( 1405 ) 
A( 1520 ) 
- 0 --- A( 1600 ) 
A( 1670 ) 
-X— A( 1690 ) 
A( 1800 ) 
■Q-- A( 1810 ) 
A( 1820 } 
A( 1830 ) 
-A--- A( 1890 } 
A( 2100 } 
-A-- A( 2110 } 


2.2 2.4 2.6 

[GeV] 



Figure 3: Fit projections for (a) rriKp and (b) mj/^p for the reduced A* model with two P+ states 
(see Table [^. The data are shown as solid (black) squares, while the solid (red) points show the 
results of the ht. The solid (red) histogram shows the background distribution. The (blue) open 
squares with the shaded histogram represent the Pc(4:450)+ state, and the shaded histogram 
topped with (purple) hlled squares represents the Pc(4380)+ state. Each A* component is also 
shown. The error bars on the points showing the ht results are due to simulation statistics. 
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Analysis and results 


We use data corresponding to 1 fb“^ of integrated luminosity acquired by the LHCb 
experiment in pp collisions at 7 TeV center-of-mass energy, and 2 fb“^ at 8 TeV. The 
LHCb detector 13 is a single-arm forward spectrometer covering the pseudorapidity range 


2 < p < 5. The detector includes a high-precision tracking system consisting of a silicon- 
strip vertex detector surrounding the pp interaction region 14 , a large-area silicon-strip 


detector located upstream of a dipole magnet with a bending power of about 4Tm, and 
three stations of silicon-strip detectors and straw drift tubes 15 placed downstream of 


the magnet. Different types of charged hadrons are distinguished using information from 
two ring-imaging Cherenkov detectors 16 . Muons are identihed by a system composed of 
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alternating layers of iron and multiwire proportional chambers 

Events are triggered by a J/ip —)■ decay, requiring two identihed muons with 

opposite charge, each with transverse momentum, pt, greater than 500 MeV. The dimuon 
system is required to form a vertex with a ht < 16, to be signihcantly displaced from 
the nearest pp interaction vertex, and to have an invariant mass within 120 MeV of the 
J/'ip mass 12 . After applying these requirements, there is a large J/'ip signal over a small 
Only candidates with dimuon invariant mass between 


background Only candidates with dimuon invariant mass between —48 MeV and 
-|-43 MeV relative to the observed J/'tp mass peak are selected, the asymmetry accounting 
for hnal-state electromagnetic radiation. 

Analysis preselection requirements are imposed prior to using a gradient Boosted 
Decision Tree, BDTG 19 , that separates the signal from backgrounds. Each track 


is required to be of good quality and multiple reconstructions of the same track are 
removed. Requirements on the individual particles include pt > 550 MeV for muons, and 
Pt > 250 MeV for hadrons. Each hadron must have an impact parameter with respect 
to the primary pp interaction vertex larger than 9, and must be positively identified in the 
particle identihcation system. The K~p system must form a vertex with < 16, as must 
the two muons from the J/ip decay. Requirements on the candidate include a vertex 

< 50 for 5 degrees of freedom, and a flight distance of greater than 1.5 mm. The vector 
from the primary vertex to the A^ vertex must align with the A^ momentum so that the 
cosine of the angle between them is larger than 0.999. Candidate pApT combinations are 
constrained to the J/A mass for subsequent use in event selection. 

The BDTG technique involves a “training” procedure using sideband data background 
and simulated signal samples. (The variables used are listed in the supplementary material.) 
We use 2 x 10® T® —)■ J/ipK~p events with J/ijj —)■ p~ that are generated uniformly in 

phase space in the LHCb acceptance, using Pythia 
tune 
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21 , and the LHCb detector simulation based on Geant4 22 


with a special LHCb parameter 
described in Ref. 
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The product of the reconstruction and trigger efficiencies within the LHCb geometric 
acceptance is about 10%. In addition, specihc backgrounds from and R® decays 
are vetoed. This is accomplished by removing combinations that when interpreted as 
J/AK^K~ fall within ±30 MeV of the R® mass or when interpreted as J/'ipK~n^ fall 
within ±30 MeV of the R® mass. This requirement effectively eliminates background 
from these sources and causes only smooth changes in the detection efficiencies across the 
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Figure 4: Invariant mass spectrum of J/ipK p combinations, with the total fit, signal and 
background components shown as solid (blue), solid (red) and dashed lines, respectively. 


decay phase space. Backgrounds from Si, decays cannot contribnte significantly to 
onr sample. We choose a relatively tight cnt on the BDTG output variable that leaves 
26 007±166 signal candidates containing 5.4% background within ±15 MeV (±2cr) of the 
J/'rpK~p mass peak, as determined by the unbinned extended likelihood ht shown in Fig. 
The combinatorial background is modeled with an exponential function and the signal 
shape is parameterized by a double-sided Hypatia function , where the signal radiative 
tail parameters are fixed to values obtained from simulation. For subsequent analysis we 
constrain the J/ijj K~p four-vectors to give the A^ invariant mass and the A^ momentum 
vector to be aligned with the measured direction from the primary to the A^ vertices 25 . 
In Fig. I^we show the “Dalitz” plot using the K~p and J/rjjp invariant masses- 


squared as independent variables. A distinct vertical band is observed in the K~p invariant 
mass distribution near 2.3 GeV^ corresponding to the 71(1520) resonance. There is also a 
distinct horizontal band near 19.5 GeV^. As we see structures in both K~p and J/tpp mass 
distributions we perform a full amplitude analysis, using the available angular variables 
in addition to the mass distributions, in order to determine the resonances present. No 
structure is seen in the J/'ipK~ invariant mass. 

We consider the two interfering processes shown in Fig. which produce two distinct 
decay sequences; A° —)■ J/^jA*, A* —)■ K~p and A° —)■ P^K~, —)■ J/i^p, with J/iIj —)■ 

pApi~ in both cases. We use the helicity formalism in which each sequential decay 
A ^ BC contributes to the amplitude a term 


KPB->.J^B,OA,oyR4mBc) = ( 1 ) 


where A is the quantum number related to the projection of the spin of the particle onto 
its momentum vector (helicity) and are complex helicity-coupling amplitudes 
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Figure 5: Invariant mass squared of K~p versus J/ij^p for candidates within ±15 MeV of the 


mass. 


describing the decay dynamics. Here 9a and (pB are the polar and azimuthal angles of B 
in the rest frame of A {6a is known as the “helicity angle” of A). The three arguments of 
Wigner’s D-matrix are Euler angles describing the rotation of the initial coordinate system 
with the 2 ;-axis along the helicity axis of A to the coordinate system with the 2 ;-axis along 
the helicity axis of B [^. We choose the convention in which the third Euler angle is 
zero. In Eq. (jl |, Wigner small-d matrix. If A has a non-negligible 

natural width, the invariant mass distribution of the B and C daughters is described by 
the complex function RA{mBc) discussed below, otherwise RA{mBc) = 1- 

Using Clebsch-Gordan coefficients, we express the helicity couplings in terms of LS 
couplings {Bb^s)i where L is the orbital angular momentum in the decay, and S is the 
total spin of A plus B\ 


njA^BC 

^Xb,Xc 


EE 


2L+1 p , 

I 


Jb Jc 

B —Ac 


5 

As — A 


c 


X 


L S 
0 Xb — Xc 


Ja 

Xb — Ac 


( 2 ) 

where the expressions in parentheses are the standard Wigner 3j-symbols. For strong decays, 
possible L values are constrained by the conservation of parity (P): Pa = PbPc{-1)^- 
Denoting J/'ijj as ip, the matrix element for the A^ -A- J/ipA* decay sequence is 


M 


E E E a,.-a,(o. Oai, 0)' 




TL A/l* \ql 


where the a:-axis, in the coordinates describing the A^ decay, is chosen to hx (pA* = 0. The 
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sum over n is due to many different resonances contributing to the amplitude. Since 
the J/-^ decay is electromagnetic, the values of AA^ = A^+ — A^- are restricted to ±1. 

There are 4 (6) independent complex couplings to £t for each resonance 

for Jyi* = I (> \)- They can be reduced to only 1 (3) free Bi s couplings to fit if only the 
lowest (the lowest two) values of L are considered. The mass rriKp, together with all decay 
angles entering Eq. (j^, 0^0, 9a*-, (j)K, 9.,p and cj)^ (denoted collectively as i7), constitute the 
six independent dimensions of the T® —)■ J/ippK~ decay phase space. 

Similarly, the matrix element for the decay chain is given by 


^ p p 


A A? 


«b<bO)E(4) 


where the angles and helicity states carry the superscript or subscript Pc to distinguish 
them from those defined for the A* decay chain. The sum over j allows for the possibility 
of contributions from more than one P^ resonance. There are 2 (3) independent helicity 

couplings Jpcj = 2 2 )’ ^ ratio of the two couplings, to 

ijj ^ 

determine from the data. 

The mass-dependent RA^irriKp) and Rp^.{mj/.,jjp) terms are given by 

Rx{m) = B^^^Jp,Po,d) BW{m\Mox,rox)B^i^{q,qo,d) (5) 


Here p is the A = H* or momentum in the A^ rest frame, and q is the momentum 
of either decay product of X in the X rest frame. The symbols po and go denote values 
of these quantities at the resonance peak {m = Mqx)- The orbital angular momentum 
between the decay products of A^ is denoted as L^o- Similarly, Lx is the orbital angular 
momentum between the decay products of X. The orbital angular momentum barrier 
factors, p^ B'i{p,po,d), involve the Blatt-Weisskopf functions 28 , and account for the 


difficulty in creating larger orbital angular momentum L, which depends on the momentum 
of the decay products p and on the size of the decaying particle, given by the d constant. 
We set d = 3.0 GeV~^ ~ 0.6 fm. The relativistic Breit-Wigner amplitude is given by 


where 


BW(m|Mox, Box) 


1 

M„x" - m" - iM„xT(m) ’ 


( 6 ) 


r(m) = r„xl-l —B,^Aq,q„,d), P) 

is the mass dependent width of the resonance. For the yl(1405) resonance, which peaks 
below the K~p threshold, we use a two-component FlattGlike parameterization (see 
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the supplementary material). The couplings for the allowed channels, Ett and Kp, are 
taken to be equal and to correspond to the nominal value of the width j^. For all 
resonances we assume minimal values of L^q and of Lx in Rxijn). For nonresonant (NR) 
terms we set BW(m) = 1 and Mqnr to the midrange mass. 

Before the matrix elements for the two decay sequences can be added coherently, the 
proton and muon helicity states in the A* decay chain must be expressed in the basis of 
helicities in the decay chain. 






b 


+ e' 


aa„ 


E^ 


1 

2 

\Pc \ 
Av 1 A'i 


% 


P 


( 8 ) 


where Op is the polar angle in the p rest frame between the boost directions from the A* 
and rest frames, and is the azimuthal angle correcting for the difference between 
the muon helicity states in the two decay chains. Note that m^p, 0^1, (l)p^, Op^, 0^, 

^b 

0^’', 0^'=, Op and can all be derived from the values of mxp and i7, and thus do not 
constitute independent dimensions in the A’^ decay phase space. (A detailed prescription 
for calculation of all the angles entering the matrix element is given in the supplementary 
material.) 

Strong interactions, which dominate A^ production at the LHC, conserve parity and 
cannot produce longitudinal A^ polarization |^. Therefore, A^o = +1/2 and —1/2 values 
are equally likely, which is reflected in Eq. (|^. If we allow the A^ polarization to vary, the 
data are consistent with a polarization of zero. Interferences between various A^ and P^^ 
resonances vanish in the integrated rates unless the resonances belong to the same decay 
chain and have the same quantum numbers. 

The matrix element given by Eq. ([^ is a 6-dimensional function of mxp and i7 and 
depends on the fit parameters, at, which represent independent helicity or LS couplings, 
and masses and widths of resonances (or Elatte parameters), AA. = AA{mKp, i7| at). After 
accounting for the selection efficiency to obtain the signal probability density function 
(PDE) an unbinned maximum likelihood fit is used to determine the amplitudes. Since the 
efhciency does not depend on at, it is needed only in the normalization integral, which is 
carried out numerically by summing |Ad(mxp, i7| at)| over the simulated events generated 
uniformly in phase space and passed through the selection. (More details are given in the 
supplementary material.) 

We use two ht algorithms, which were independently coded and which differ in the 
approach used for background subtraction. In the hrst approach, which we refer to as 
cEit, the signal region is defined as ±2 a around the A^ mass peak. The total PDE used 
in the ht to the candidates in the signal region, V{mKp, f2| at), includes a background 
component with normalization hxed to be 5.4% of the total. The background PDF is 
found to factorize into hve two-dimensional functions of rriKp and of each independent 
angle, which are estimated using sidebands extending from 5.0 a to 13.5 a on both sides of 
the peak. 

In the complementary approach, called sFit, no explicit background parameterization is 
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needed. The PDF consists of only the signal component, with the background subtracted 
using the sPlot technique 31 applied to the log-likelihood sum. All candidates shown in 
Fig. I^are included in the sum with weights, Wi, dependent on mj/^Kp- The weights are 
set according to the signal and the background probabilities determined by the hts to the 
mj/^pK distributions, similar to the ht displayed in Fig. but performed in 32 different 
bins of the two-dimensional plane of cos Oy^o and cos 6j/^ to account for correlations with 
the mass shapes of the signal and background components. This quasi-log-likelihood 
sum is scaled by a constant factor, sw = fo account for the effect of the 

background subtraction on the statistical uncertainty. (More details on the cFit and sFit 
procedures are given in the supplementary material.) 

In each approach, we minimize —21n£(at) = —2svj/ITj lnP(m/i'p., i7j| at), which 
gives the estimated values of the ht parameters, at min, together with their covariance 
matrix [Wi = 1 in cFit). The difference of —21n£(atmin) between different amplitude 
models, A(—21n£), allows their discrimination. For two models representing separate 
hypotheses, e.g. when discriminating between different values assigned to a Pf state, 
the assumption of a distribution with one degree of freedom for A(—21n£) under the 
disfavored hypothesis allows the calculation of a lower limit on the signihcance of its 
rejection, i.e. the p-value 32 . Therefore, it is convenient to express A(—21n£) values as 


n^, where Uo- corresponds to the number of standard deviations in the normal distribution 
with the same p-value. For nested hypotheses, e.g. when discriminating between models 
without and with P+ states, overestimates the p-value by a modest amount. Simulations 
are used to obtain better estimates of the signihcance of the P+ states. 

Since the isospin of both the and the J/^jj particles are zero, we expect that the 
dominant contributions in the K~p system are A* states, which would be produced via a 
AJ = 0 process. It is also possible that E* resonances contribute, but these would have 
AJ = 1. By analogy with kaon decays the AJ = 0 process should be dominant 33 . The 
list of A* states considered is shown in Table [U 

Our strategy is to hrst try to ht the data with a model that can describe the mass and 
angular distributions including only A* resonances, allowing all possible known states and 
decay amplitudes. We call this the “extended” model. It has 146 free parameters from the 
helicity couplings alone. The masses and widths of the A* states are hxed to their PDG 
values, since allowing them to hoat prevents the ht from converging. Variations in these 
parameters are considered in the systematic uncertainties. 

The cFit results without any P+ component are shown in Fig. While the rrixp 
distribution is reasonably well htted, the peaking structure in mj/^ijp is not reproduced. 
The same result is found using sFit. The speculative addition of E* resonances to the 
states decaying to K~p does not change this conclusion. 

We will demonstrate that introducing two —)■ J/tpp resonances leads to a satisfactory 

description of the data. When determining parameters of the states, we use a more 
restrictive model of the K~p states (hereafter referred to as the “reduced” model) that 
includes only the A* resonances that are well motivated, and has fewer than half the number 
of free parameters. As the minimal L(Jo for the spin 9/2 A(2350) equals Ja* — Ja^ ~ Jj/ip = 3, 
it is extremely unlikely that this state can be produced so close to the phase space limit. 
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Table 1: The A* resonances used in the different fits. Parameters are taken from the PDG 
We take 5/2“ for the of the 71(2585). The number of LS couplings is also listed for both 
the “reduced” and “extended” models. To fix overall phase and magnitude conventions, which 
otherwise are arbitrary, we set i = (1,0) for 71(1520). A zero entry means the state is excluded 
from the fit. 


State 

JP 

Mo (MeV) 

To (MeV) 

# Reduced 

# Extended 

71(1405) 

1/2- 

1405.ll|;|j 

50.5 ±2.0 

3 

4 

71(1520) 

3/2“ 

1519.5 ± 1.0 

15.6 ± 1.0 

5 

6 

71(1600) 

1/2+ 

1600 

150 

3 

4 

71(1670) 

1/2- 

1670 

35 

3 

4 

71(1690) 

3/2“ 

1690 

60 

5 

6 

71(1800) 

1/2- 

1800 

300 

4 

4 

71(1810) 

1/2+ 

1810 

150 

3 

4 

71(1820) 

5/2+ 

1820 

80 

1 

6 

71(1830) 

5/2“ 

1830 

95 

1 

6 

71(1890) 

3/2+ 

1890 

100 

3 

6 

71(2100) 

7/2“ 

2100 

200 

1 

6 

71(2110) 

5/2+ 

2110 

200 

1 

6 

71(2350) 

9/2+ 

2350 

150 

0 

6 

71(2585) 

? 

^2585 

200 

0 
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Figure 6: Results for (a) rriKp and (b) rrij/^p for the extended 7l* model fit without states. 
The data are shown as (black) squares with error bars, while the (red) circles show the results of 
the fit. The error bars on the points showing the fit results are due to simulation statistics. 
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In fact L = 3 is the highest orbital angular momentum observed, with a very small rate, 
in decays of B mesons with much larger phase space available {Q = 2366 MeV, while 
here Q = 173 MeV), and without additional suppression from the spin counting factors 
present in 71(2350) production (all three Jyi*, J^o and Jj/^ vectors have to line up in the 

same direction to produce the minimal value). Therefore, we eliminate it from the 

0 

reduced A* model. We also eliminate the 7l(2585) state, which peaks beyond the kinematic 
limit and has unknown spin. The other resonances are kept but high Ldo amplitudes are 
removed; only the lowest values are kept for the high mass resonances, with a smaller 
reduction for the lighter ones. The number of LS amplitudes used for each resonance is 
listed in Table With this model we reduce the number of parameters needed to describe 
the A* decays from 146 to 64. For the different combinations of resonances that we 
try, there are up to 20 additional free parameters. Using the extended model including 
one resonant improves the ht quality, but it is still unacceptable (see supplementary 
material). We hnd acceptable hts with two states. We use the reduced A* model for 
the central values of our results. The differences in htted quantities with the extended 
model are included in the systematic uncertainties. 

The best ht combination hnds two states with values of 3/2“ and 5/2’'', for the 
lower and higher mass states, respectively. The —2 ln£ values differ by only 1 unit between 
the best ht and the parity reversed combination (3/2’'', 5/2“). Other combinations are 
less likely, although the (5/2’'', 3/2“) pair changes —21n£ by only 2.3^ units and therefore 
cannot be ruled out. All combinations 1/2^ through 7/2^ were tested, and all others 
are disfavored by changes of more than 5^ in the —21n£ values. The cFit results for the 
(3/2“, 5/2+) ht are shown in Fig. Both distributions of rriKp and mj/^p are reproduced. 
The lower mass 3/2“ state has mass 4380±8 MeV and width 205±18 MeV, while the 5/2+ 
state has a mass of 4449.8±1.7 MeV and width 39±5 MeV; these errors are statistical only, 
systematic uncertainties are discussed later. The mass resolution is approximately 2.5 MeV 
and does not affect the width determinations. The sFit approach gives comparable results. 
The angular distributions are reasonably well reproduced, as shown in Fig. and the 
comparison with the data in rriKp intervals is also satisfactory as can be seen in Fig. 
Interference effects between the two P+ states are particularly evident in Figj^d), where 
there is a large destructive contribution (not explicitly shown in the hgure) to the total 
rate. (A ht fraction comparison between cFit and sFit is given in the supplementary 
material.) The addition of further P+ states does not signihcantly improve the ht. 

Adding a single 5/2+ P+ state to the ht with only A* states reduces —21n£ by 14.7^ 
using the extended model and adding a second lower mass 3/2“ P+ state results in a further 
reduction of 11.6^. The combined reduction of —21n£ by the two states taken together is 
18.7^. Since taking a/A 2 In C overestimates signihcances, we perform simulations to obtain 
more accurate evaluations. We generate pseudoexperiments using the null hypotheses 
having amplitude parameters determined by the hts to the data with no or one P+ state. 
We ht each pseudoexperiment with the null hypothesis and with P+ states added to the 
model. The —2 In C distribntions obtained from many psendoexperiments are consistent 
with distributions with the number of degrees of freedom approximately equal to twice 
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Figure 7: Various decay angular distributions for the fit with two states. The data are shown 
as (black) squares, while the (red) circles show the results of the fit. Each fit component is also 
shown. The angles are defined in the text. 
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Figure 8: mj/^p in various intervals of rriKp for the fit with two states: (a) rriKp < 1-55 GeV, 
(b) 1.55 < rriKp < 1-70 GeV, (c) 1.70 < rriKp < 2.00 GeV, and (d) rriKp > 2.00 GeV. The data 
are shown as (black) squares with error bars, while the (red) circles show the results of the fit. 
The blue and purple histograms show the two states. See Fig. Mfor the legend. 
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the number of extra parameters in the ht. Comparing these distributions with the A2 ln£ 
values from the hts to the data, p-values can be calculated. These studies show reduction 
of the signihcances relative to VA21n£ by about 20%, giving overall signihcances of 9 a 
and 12 a, for the lower and higher mass states, respectively. The combined signihcance 
of two states is 15 a. Use of the extended model to evaluate the significance includes 
the effect of systematic uncertainties due to the possible presence of additional A* states 
or higher L amplitudes. 

Systematic uncertainties are evaluated for the masses, widths and £t fractions of the 
states, and for the ht fractions of the two lightest and most signihcant A* states. 
Additional sources of modeling uncertainty that we have not considered may affect the 
ht fractions of the heavier A* states. The sources of systematic uncertainties are listed 
in Table They include diherences between the results of the extended versus reduced 
model, varying the A* masses and widths, uncertainties in the identihcation requirements 
for the proton, and restricting its momentum, inclusion of a nonresonant amplitude in 
the ht, use of separate higher and lower A^ mass sidebands, alternate hts, varying 
the Blatt-Weisskopf barrier factor, d, between 1.5 and 4.5 GeV“^, changing the angular 
momentum L used in Eq. ([^ by one or two units, and accounting for potential mismodeling 
of the efficiencies. For the /1(1405) ht fraction we also added an uncertainty for the Flatte 
couplings, determined by both halving and doubling their ratio, and taking the maximum 
deviation as the uncertainty. 

The stability of the results is cross-checked by comparing the data recorded in 2011/2012, 
with the LHCb dipole magnet polarity in up/down conhgurations, A^/A^ decays, and A^ 
produced with low/high values of px- Extended model hts without including P^ states 
were tried with the addition of two high mass A* resonances of freely varied mass and 
width, or four nonresonant components up to spin 3/2; these do not explain the data. The 
htters were tested on simulated pseudoexperiments and no biases were found. In addition, 
selection requirements are varied, and the vetoes of and are removed and explicit 
models of those backgrounds added to the ht; all give consistent results. 

Further evidence for the resonant character of the higher mass, narrower, P^ state is 
obtained by viewing the evolution of the complex amplitude in the Argand diagram (I^ . 
In the amplitude hts discussed above, the Pc(4:450) is represented by a Breit-Wigner 
amplitude, where the magnitude and phase vary with mj/^p according to an approximately 
circular trajectory in the (ReA^'=, ImA^'=) plane, where A^‘^ is the mj/^p dependent 
part of the Pc(4450)+ amplitude. We perform an additional ht to the data using the 
reduced A* model, in which we represent the Pc(4450)+ amplitude as the combination of 
independent complex amplitudes at six equidistant points in the range iFo = 39MeV 
around Mq = 4449.8 MeV as determined in the default ht. Real and imaginary parts of 
the amplitude are interpolated in mass between the htted points. The resulting Argand 
diagram, shown in Fig. [^a), is consistent with a rapid counter-clockwise change of the 
Pc(4450)+ phase when its magnitude reaches the maximum, a behavior characteristic of 
a resonance. A similar study for the wider state is shown in Fig. [^b); although the ht 
does show a large phase change, the amplitude values are sensitive to the details of the A* 
model and so this latter study is not conclusive. 
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Different binding mechanisms of pentaquark states are possible. Tight-binding was 
envisioned originally (3p,35 . A possible explanation is heavy-light diquarks (3^. Examples 


of other mechanisms include a diquark-diquark-antiquark model 37,38 , a diquark-triquark 
model 39 , and a coupled channel model 40 . Weakly bound “molecules” of a baryon plus 


41 


a meson have been also discussed 

Models involving thresholds or “cusps” 


candidates via nonresonant scattering mechanisms 42 44 


have been invoked to explain some exotic meson 
There are certain obvious 


difficulties with the use of this approach to explain our results. The closest threshold to 
the high mass state is at 4457.1±0.3 MeV resulting from a Ac{2595)~^D^ combination, 
which is somewhat higher than the peak mass value and would produce a structure with 
quantum numbers = 1/2’'' which are disfavored by our data. There is no threshold 
close to the lower mass state. 

In conclusion, we have presented a full amplitude ht to the —)■ J/'iljK~p decay. 
We observe signihcant A* production recoiling against the J/V’ with the lowest mass 
contributions, the 71(1405) and 71(1520) states having ht fractions of (15 ± 1 ± 6)% and 
(19±1±4)%, respectively. The data cannot be satisfactorily described without including two 
Breit-Wigner shaped resonances in the J/'ifjp invariant mass distribution. The signihcances 
of the lower mass and higher mass states are 9 and 12 standard deviations, respectively. 


Table 2: Summary of systematic uncertainties on masses, widths and fit fractions, and 7l* fit 
fractions. A fit fraction is the ratio of the phase space integrals of the matrix element squared 
for a single resonance and for the total amplitude. The terms “low” and “high” correspond to 
the lower and higher mass states. The sFit/cFit difference is listed as a cross-check and not 
included as an uncertainty. 


Source 

Mo (MeV) 
low high 

Bo 

low 

(MeV) 

high 

low 

Fit fractions (%) 
high A(1405) 7l(1520) 

Extended vs. reduced 

21 

0.2 

54 

10 

3.14 

0.32 

1.37 

0.15 

A* masses & widths 

7 

0.7 

20 

4 

0.58 

0.37 

2.49 

2.45 

Proton ID 

2 

0.3 

1 

2 

0.27 

0.14 

0.20 

0.05 

10 < pj, < 100 GeV 

0 

1.2 

1 

1 

0.09 

0.03 

0.31 

0.01 

Nonresonant 

3 

0.3 

34 

2 

2.35 

0.13 

3.28 

0.39 

Separate sidebands 

0 

0 

5 

0 

0.24 

0.14 

0.02 

0.03 

JP (3/2+, 5/2-) or (5/2+, 3/2“) 

10 

1.2 

34 

10 

0.76 

0.44 



d = 1.5 - 4.5 GeV-i 

9 

0.6 

19 

3 

0.29 

0.42 

0.36 

1.91 

Al P+ (low/high) A“ 

6 

0.7 

4 

8 

0.37 

0.16 



Lp^ (low/high) — J/'tpp 

4 

0.4 

31 

7 

0.63 

0.37 



if, Al ^ 

11 

0.3 

20 

2 

0.81 

0.53 

3.34 

2.31 

Efficiencies 

1 

0.4 

4 

0 

0.13 

0.02 

0.26 

0.23 

Change 7l(1405) coupling 

0 

0 

0 

0 

0 

0 

1.90 

0 

Overall 

29 

2.5 

86 

19 

4.21 

1.05 

5.82 

3.89 

sFit/cFit cross check 

5 

1.0 

11 

3 

0.46 

0.01 

0.45 

0.13 


14 




















Figure 9: Fitted values of the real and imaginary parts of the amplitudes for the baseline (3/2“, 
5/2''') fit for a) the Pc(4450)''' state and b) the Pc(4380)'’“ state, each divided into six bins 

of equal width between —Fq and +ro shown in the Argand diagrams as connected points with 
error bars {mjj^p increases counterclockwise). The solid (red) curves are the predictions from 
the Breit-Wigner formula for the same mass ranges with Mq (Fq) of 4450 (39) MeV and 4380 
(205) MeV, respectively, with the phases and magnitudes at the resonance masses set to the 
average values between the two points around Mq. The phase convention sets B^i = (1,0) for 
yl(1520). Systematic uncertainties are not included. 


These structures cannot be accounted for by reflections from J/ipA* resonances or other 
known sources. Interpreted as resonant states they must have minimal quark content of 
ccuud, and would therefore be called charmonium-pentaquark states. The lighter state 
Fc( 4380)’'' has a mass of 4380 ± 8 ± 29 MeV and a width of 205 ± 18 ± 86 MeV, while the 
heavier state Pc(4450)+ has a mass of 4449.8 ± 1.7 ± 2.5 MeV and a width of 39 ± 5 ± 19 
MeV. A model-independent representation of the Pc(4450)’'' contribution in the ht shows 
a phase change in amplitude consistent with that of a resonance. The parities of the two 
states are opposite with the preferred spins being 3/2 for one state and 5/2 for the other. 
The higher mass state has a fit fraction of (4.1 ± 0.5 ± 1.1)%, and the lower mass state of 
(8.4 ± 0.7 ± 4.2)%, of the total —)■ J/'iljK~p sample. 
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1 Variables used in the BDTG 


Muon identification uses information from several parts of the detector, including the 
RICH detectors, the calorimeters and the muon system. Likelihoods are formed for the 
muon and pion hypotheses. The difference in the logarithms of the likelihoods, DLL(p, —tt), 
is used to distinguish between the two 16 . The smaller valne of the two discriminants 


DLL(/r’'' — TT^) and DLL(/r“ — tt~) is used as one of the BDTG variables. 

The next set of variables uses the kaon and proton tracks. The Xip is dehned as the 
difference in of the primary vertex reconstructed with and without the considered track. 
The smaller Xip of the K~ and p is used in the BDTG. The scalar px sum of the K~ and 
p is another variable. 

The last set of variables uses the candidate. The cosine of the angle between a 
vector from the primary vertex to the vertex and the A^ momentnm vector is one 
input variable. In addition the Xjp, the flight distance, the px and the vertex of the A^ 
candidate are used. 


2 Additional fit results 


2.1 Reduced model fit projections for 


The Dalitz plots for the other two possible projections are shown in Fig. 10 
obvious resonance structure in the J/'ipK 


There is no 


mass-squared distribution. 

Onr ht describes well the distribution as shown by viewing the projections of 


the reduced model fit. They are shown for different slices of rriKp in Fig. 11 
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Figure 10: (a) Invariant mass squared of J/'ipK versus J/ipp and (b) of J/ipK versus K p for 
candidates within ±15 MeV of the A!/ mass. 
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Figure 11: Projections onto rrij/^x iii various intervals of rriKp for the reduced model fit (cFit) with 
two states of equal to 3/2“ and 5/2+: (a) rriKp < 1-55 GeV, (b) 1.55 < niKp < 1-70 GeV, 
(c) 1.70 < rriKp < 2.00 GeV, (d) rriKp > 2.00 GeV, and (e) all niKp- The data are shown as 
(black) squares with error bars, while the (red) circles show the results of the ht. The individual 
resonances are given in the legend. 


2.2 Reduced model angular fits with two states for 

m{K~p) > 2 GeV 

In Fig. [^we show the result of the reduced model ht to the angular distributions for 
m{K~p) > 2 GeV. The data is well described by the hts. 
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Figure 12: Various decay angular distributions for the fit with two states for m{K~p) > 2 GeV. 
The data are shown as (black) squares, while the (red) circles show the results of the fit. Each fit 
component is also shown. The angles are defined in the text. 
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Figure 13: Results of the fit with one = 5/2^ candidate, (a) Projection of the invariant 
mass of K~p combinations from —>• J/'il)K~p candidates. The data are shown as (black) 

squares with error bars, while the (red) circles show the results of the fit; (b) the corresponding 
J/il^p mass projection. The (blue) shaded plot shows the projection, the other curves represent 
individual A* states. 


2.3 Extended model fit with one 


In the fits with one P/ amplitnde, we test valnes of 1/2^, 3/2^ and 5/2^. The mass 
and width of the putative P+ state are allowed to vary. There are a total of 146 free 
parameters for the A* states to which we add either three complex couplings for 1/2^ or 
four for higher spins. The best ht is with a 5/2’'' state, which improves —21n£ by 215. 


Figure 13 shows the projections for this ht. While the iriKp projection is well described, 
clear discrepancies in mj/^p remain visible. 


2.4 Results of extended model fit with two Pc^ states 

For completeness we include here the results of the extended model ht with two 
states using cFit. We hnd acceptable hts for several combinations. For a lower mass 
= 3/2“ state and a higher mass 5/2+ state, the masses (widths) are 4358.9±6.6 MeV 
(151.1±13.7 MeV), and 4450.1±1.7 MeV (48.6±4.0 MeV), respectively. The uncertainties 


are statistical only. The results for this two P+ ht are shown in Fig. 14 Both the rriKp 
distribution and the peaking structure in mj/^p are reproduced. 


3 Fit fraction comparison between cFit and sFit 

The ht fraction for a given resonance is a ratio of the phase space integrals of the matrix 
element squared calculated for the resonance amplitude taken alone and for the total 
matrix element summing over all contributions. The ht fractions are listed in Table The 
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Figure 14: Results from cFit for (a) rriKp and (b) rrij/^p for the extended model with two 
states. The data are shown (black) squares with error bars, while the (red) circles show the 
results of the fit. Each A* component is also shown. The (blue) open squares and (purple) solid 
squares show the two states. 


states have well determined fit fractions. There is good agreement between cFit and 
sFit. Note that the results for the yl(1405) resonance are based on our use of a particular 
Flatte amplitude model. 

Table 3: Fit fractions of the different components from cFit and sFit for the default (3/2“, 5/2~^) 
model. Uncertainties are statistical only. 


Particle Fit fraction (%) cFit Fit fraction (%) sFit 

Pc (4380) + 

8.42 ±0.68 

7.96 ±0.67 

Pc (4450) + 

4.09 ±0.48 

4.10 ±0.45 

T(1405) 

14.64 ±0.72 

14.19 ±0.67 

T(1520) 

18.93 ±0.52 

19.06 ±0.47 

T(1600) 

23.50 ± 1.48 

24.42 ± 1.36 

T(1670) 

1.47 ±0.49 

1.53 ±0.50 

T(1690) 

8.66 ±0.90 

8.60 ±0.85 

T(1800) 

18.21 ±2.27 

16.97 ±2.20 

21(1810) 

17.88 ±2.11 

17.29 ±1.85 

21(1820) 

2.32 ±0.69 

2.32 ±0.65 

21(1830) 

1.76 ±0.58 

2.00 ±0.53 

21(1890) 

3.96 ±0.43 

3.97 ±0.38 

21(2100) 

1.65 ±0.29 

1.94 ±0.28 

21(2110) 

1.62 ±0.32 

1.44 ±0.28 
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4 Details of the matrix element for the decay ampli¬ 
tude 

The matrix element for —>■ 'iljK~p, 'll; —)■ decay^ must allow for various conven¬ 
tional A* —)■ K~p resonances and exotic pentaquark states —)■ ipp that could interfere 

with each other. 

We use the helicity formalism to write down the matrix element. To make the derivation 
of the matrix element easier to comprehend we start with a brief outline of this formalism 
and our notation. Then we discuss the application to the —)■ A* —)■ K~p, 

'ip —)■ pApi~ decay sequence, called hereafter the A* decay chain matrix element. Next we 
discuss construction of the A^ —)■ Pp'K~, ipp, tp —)■ decay sequence, called 

hereafter the Pc decay chain matrix element, which can be coherently added to that for the 
A* decay chain. We also discuss a possible reduction of the number of helicity couplings 
to be determined from the data using their relationships to the LS couplings. 

4.1 Helicity formalism and notation 

For each two-body decay A ^ B C, & coordinate system is set up in the rest frame of A, 
with 5 bein^the direction of quantization for its spin. We denote this coordinate system 
as where the superscript “{T}” means “in the rest frame of A”, while 

the subscript “0” means the initial coordinates. For the first particle in the decay chain 
(4°), the choice of these coordinates is arbitrary]^ However, once defined, these coordinates 
must be used consistently between all decay sequences described by the matrix element. 
For subsequent decays, e.g. B ^ D E, the choice of these coordinates is already fixed by 
the transformation from the A to the B rest frames, as discussed below. Helicity is dehned 
as the projection of the spin of the particle onto the direction of its momentum. When 
the z axis coincides with the particle momentum, we denote its spin projection onto it 
{i.e. the rriz quantum number) as A. To use the helicity formalism, the initial coordinate 
system must be rotated to align the 2 ; axis with the direction of the momentum of one of 
the daughter particles, e.g. the B. A generalized rotation operator can be formulated in 
three-dimensional space, that uses Euler angles. Applying this operator results 

in a sequence of rotations: first by the angle a about the Zq axis, followed by the angle (3 
about the rotated yi axis and then hnally by the angle 7 about the rotated Z 2 axis. We use 
a subscript denoting the axes, to specify the rotations which have been already performed 
on the coordinates. The spin eigenstates of particle A, in the {xq ,yfj^, Zq"^) 

coordinate system can be expressed in the basis of its spin eigenstates, \JA,m'^), in the 

^We denote J/p as p for efficiency of the notation. 

^The “hat” symbol denotes a unit vector in a given direction. 

“^When designing an analysis to be sensitive (or insensitive) to a particular case of polarization, the 
choice is not arbitrary, but this does not change the fact that one can quantize the spin along any 
well-defined direction. The polarization may be different for different choices. 
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( 9 ) 


rotated 


coordinate system with the help of Wigner’s D—matrices 


\JA,'mA) = '^D 


Ja 




where 


= {J,m\n{a,P,j)\J,mr = e*”^“ d 


J 

m,m‘ 


m e'"'-', 


( 10 ) 


and where the small-d Wigner matrix contains known functions of (3 that depend on 
To achieve the rotation of the original Zq^ axis onto the B momentum 
it is sufficient to rotate by a = (3 = 0^ ^ where ■, 0^ are the azimuthal and 

polar angles of the B momentum vector in the original coordinates i.e. 

This is depicted in Fig. for the case when the quantization axis for the spin of A is its 
momentum in some other reference frame. Since the third rotation is not necessary, we 
set 7 = 00 The angle 6^"^ is usually called “the A helicity angle”, thus to simplify the 
notation we will denote it as 6a- For compact notation, we will also denote as 0 b. 
These angles can be determined fromj^ 


0 B = atan 2 

= atan 2 j 

= atan 2 x Xg^^^) ■ pj^\ Xg^^^^ ■ pj^^ j , ( 11 ) 

COsd A = ■ Pb"^^ ■ ( 12 ) 


Angular momentum conservation requires + m0 = As — Ac (since p^^ points 

in the opposite direction to Z 3 ^\ rn'(j = —Ac). Each two-body decay adds a multiplicative 
term to the matrix element 


njA^BC 


D 


Ja 

BT-A, ^B — ^C 


(<pB,OA,oy. 


(13) 


The helicity couplings are complex constants. Their products from subsequent 

decays are to be determined by the fit to the data (they represent the decay dynamics). If 
the decay is strong or electromagnetic, it conserves parity which reduces the number of 
independent helicity couplings via the relation 


njA-A-BC _ P. P P ( -\\J b+Jc—Ja ^ A-aB C 

^ — Xb,—Xc -^A J^B \ J-J ’ 


(14) 


where P stands for the intrinsic parity of a particle. 

®An alternate convention is to set 7 = —a. The two conventions lead to equivalent formulae. 

®The function atan2(x,y) is the tan~^{y/x) function with two arguments. The purpose of using 
two arguments instead of one is to gather information on the signs of the inputs in order to return the 
appropriate quadrant of the computed angle. 
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After multiplying terms given by Eq. (13) for all decays in the decay sequence, they 
must be summed up coherently over the helicity states of intermediate particles, and 
incoherently over the helicity states of the initial and final-state particles. Possible helicity 
values of B and C particles are constrained by |Ab| < Jb, |Ac| < Jc and |As — Ac| < Ja- 
When dealing with the subsequent decay of the daughter, B ^ D E, four-vectors of 
all particles must be first Lorentz boosted to the rest frame of B, along the i.e. 
direction (this is the axis in the rest frame of A after the Euler rotations; we use the 
subscript “3” for the number of rotations performed on the coordinates, because of the 
three Euler angles, however, since we use the 7 = 0 convention these coordinates are the 
same as after the first two rotations). This is visualized in Fig. 15, with B ^ D E particle 
labels replaced hy A ^ B C labels. This transformation does not change vectors that 
are perpendicular to the boost direction. The transformed coordinates become the initial 


Rest frame of A 


Helicity frame of A 



Figure 15: Coordinate axes for the spin quantization of particle A (bottom part), chosen to be 
the helicity frame of A (Aolllbl in the rest frame of its mother particle or in the laboratory frame), 
together with the polar {9^) and azimuthal angles of the momentum of its daughter B 

in the A rest frame (top part). Notice that the directions of these coordinate axes, denoted as 
, and Zq^, do not change when boosting from the helicity frame of A to its rest frame. 
After the Euler rotation TZ{a = = 9^^= 0) (see the text), the rotated z axis, , is 

aligned with the B momentum; thus the rotated coordinates become the helicity frame of H. If 
B has a sequential decay, then the same boost-rotation process is repeated to define the helicity 
frame for its daughters. 
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coordinate system quantizing the spin of B in its rest frame, 


X, 


{5} 


= X 


{4 


. {B} . (4 

Vo =% ) 


(4 


= 


(4 


( 15 ) 


The processes of rotation and subsequent boosting can be repeated until the hnal-state 
particles are reached. In practice, there are two equivalent ways to determine the 


14 


direction. Using Eq. (15) we can set it to the direction of the B momentum in the A rest 
frame 

(16) 

Alternatively, we can make use of the fact that B and C are back-to-back in the rest frame 
of A, = —p^^- Since the momentum of C is antiparallel to the boost direction from 
the A to B rest frames, the C momentum in the B rest frame will be different, but it will 
still be antiparallel to this boost direction 


;|{4 _ 


- {B} 
-Pc ■ 


(17) 


To determine x^j^^ from Eq. (15), we need to hnd x^"^^ ■ After the hrst rotation by 0s 

about Zfj^\ the x}^ axis is along the component of p^ which is perpendicular to the 
. {A} 

Zn axis 


14 


''BXzo 


= {p 


B 


-14 
= Pb 


' Xz “ Pb \bb 

J-Zq 


- (pT), 




{P 


14 

B 


{45 i (4 


A {4 _ ^5 (4 

%Jb ^ — LL 


B±zo 


After the second rotation by 9a about y-p , z. 


44 

’'BXzo 

^(4 

’'BXzo 

(4 *14 _ -(4 


(18) 


= zz ' — 


antiparallel to the component of the Zq 
i.e. Pb^^. Thus 


(4 


Pb^\ and £ 2 ^^^ = £ 3 ^^^ 


IS 


vector that is perpendicular to the new z axis 


^zo±B — 


(45 ^ (4 

0 - 4 


X, 


(4 


(4 


= xz ^ = 






^(4 


(19) 


Then we obtain = Zg^^^ x 

If C also decays, C ^ F G, then the coordinates for the quantization of C spin in the 
C rest frame are dehned by 


X, 


(4 _ (4 _ ^ {4 _ (4 

) — '^3 — Pc ~ Pb ^ 

{4 _ £ (4 _ _ ^ (4 _ , ^ (4 

0 “ -^3 “ ^zqXB — 

{4 _ (4 ^ - (4 


Vo = zp ‘ XX, 


( 20 ) 

( 21 ) 

( 22 ) 
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i.e. the ^ axis is reflected compared to the system used for the decay of particle B (it must 
point in the direction of C momentnm in the A rest frame), bnt the x axis is kept the 


same, since we chose particle B for the rotation used in Eq. (13). 


4.2 Matrix element for the A* decay chain 

We first discuss the part of the matrix element describing conventional ^ ^ Kp 

decays {i.e. A* decay chain), where denotes various possible excitations of the 4, e.g. 
4(1520). For simplicity we often refer to as 4*, unless we label an n-dependent quantity. 
The weak decay 4° —>■ 4^-^ is described by 


'IJ b 

\a* 5 




X 

b 




(23) 


where B 




are resonance {i.e. n) dependent helicity couplings to be determined by a 


fit to the data. There are 4 different complex valnes of these conplings to be determined 
for each resonance with spin = \) and 6 values for higher spins. The couplings are 
complex parameters; thus each independent coupling contributes 2 free parameters (taken 
to be real and imaginary parts) to the fit. Since the ify and A* are intermediate particles 
in the decay chain, the matrix element terms for different valnes of and A^i* mnst be 
added coherently. 

The choice of the direction for the 4° spin qnantization is arbitrary. We choose 

the 4° momentum in the lab frame to define the direction, giving its spin projection 
onto this axis the meaning of the 4° helicity (Ayig). In the 4° rest frame, this direction is 
defined by the direction of the boost from the lab frame (Eq. 


^ {lab} 
PAi 


Zn = 


(24) 


as depicted in Fig. 
as 


16 


With this choice, 0^0 is the 4° helicity angle and can be calculated 


^ « {lab} 

COS 0^0 =p^o 


Pi-"'- (25) 

Longitudinal polarization of the 4° via strong production mechanisms is forbidden due to 
parity conservation in strong interactions, cansing Ayjo = and — | to be equally likely. 

laO} 

Terms with different A^o values must be added incoherently. The choice of Xq ^ direction 
in the 4° rest frame is also arbitrary. We use the 4° —)■ 4*'0 decay plane in the lab frame 
to define it, which makes the (j)A* angle zero by definition. 

The strong decay —)■ Kp is described by a term 


-OA.f,A,('('*■.9.1-.0)- -RAjlm/tp). 


(26) 


Since the K meson is spinless, the resonance-dependent helicity coupling depends 


only on proton helicity, Ap = ±|. 
couplings are related 


H'VA” = -P4- 


As strong decays conserve parity, the two helicity 

^27) 
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Figure 16: Definition of the decay angles in the A* decay chain. 


where Pa* is the parity of A^. Since the overall magnitude and phase of can 

■*"2 

be absorbed into a redehnition of the couplings, we set 


= (1, 0) and 


2 ^0), where the values in parentheses give the real and imaginary 

“ 2 

parts of the couplings. 

The angles and 9a* are the azimuthal and polar angles of the kaon in the yl* rest 
frame (see Fig. 16). The ^ direction is defined by the boost direction from the A^ rest 
frame, which coincides with the ^ direction in this frame (Eq. (17)). This leads to 


(28) 


/J /V {A*} ^ {^*} 

cosOa* = -p^ -Pk , 

with both vectors in the A* rest frame. As explained in Sec. 
dehned by the choice of coordinates in the rest frame discussed above. Following 
Eq. (p^ and (pij), we have 


4.1 


the Xr!^ ^ direction i 


IS 




{lab} . {ag}. (ylgl 


f V*bl _ ^ - l^au/ . ^ V*6/\ - v^bf 

^zo±A* — Pao ^Pao Fa* ) Pa* 


X, 


D*} 




= Xn " = 


^ 2 o±yl* 

^zo±A* 


(29) 


The azimuthal angle of the K can now be determined in the A* rest frame from (Eq. 0 ) 


<t>K = atim2 ( - (pi'"* X • pP, ■ p};' 


(30) 


The term RA^irriKp) describes the resonance that appears in the invariant mass 
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distribution of the kaon-proton system, 


Ra^ {rriKp) = B' ^ (p, po, d) 


p 


rAt 


"^9 


M 




BW(mKp|Mo^*, rf) (g, go, d) 


"At 


M 


At 


Here, p is the A* momentum in the rest frame (p = |p^f^|)- Similarly, g is the K 

(A *1 

momentnm in the A* rest frame (g = \px ^D- The symbols po and go denote valnes of 
these qnantities at the resonance peak {niKp = M^). The orbital angnlar momentum 
between the A* and ip particles in the decay is denoted as Similarly, Lyi* is the 
orbital angular momentum between the p and K~ in the decay. For the orbital angular 


(31) 


momentum barrier factors, i?^(p, po, d), we use the Blatt-Weisskopf functions 45 
Bo{p,Po,d) =1, 

B[{p,po,d) = 


'l + (po d )2 


1 + (p d)2 ’ 



' 9 + 3(po dy + (po dY 
9 + 3(p dY + (p dY ’ 


(32) 


'225 + 45(po dy + 6(po dy + (po d)^ 
225 + 45(p dy + 6(p dy + (p dy 


H1025 + 1575(po dy + 135(po dy + 10(po dy + (po dy 
11025 + 1575(p d)2 + 135(p d)^ + 10(p dy + {p dy ’ 

1 893025 + 99225(po dy + 6300(po dy + 315(po dy + 15(po dy + (po dy^ 
893025 + 99225(p d)2 + 6300(p d)4 + 315(p d)6 + 15(p d)® + (p ’ 


to account for the difficulty in creating the orbital angular momentum L, which depends 
on the momentnm of the decay prodncts p (in the rest frame of the decaying particle) 
and on the size of the decaying particle given by the constant d. We set d = 3.0 GeV“^ 
~0.6 fm. The relativistic Breit-Wigner amplitude is given by 


BW(m|Mo,ro) 


1 

m 2 - m2 - iMoT{m) ’ 


(33) 


where 

r(m) = Fo (g, go, df . (34) 

\qoJ rn 

In the case of the 71(1405) resonance, which peaks below the K~p threshold, we 
use a two-component width equivalent to the Flatte parameterization. We add to the 
above width in the K~p channel, a width for its decay to the dominant channel, 

F(m) = T{m)K-p + I'{rn)z]n, where g in the second term and go in both terms are calculated 
assuming the decay to . Assuming that both channels are dynamically equally likely 
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and differ only by the phase space factors we set Fq to the total width of yl(1405) in both 
terms. 

Angular momentum conservation limits L 4 * to J/i* ± which is then uniquely dehned 
by parity conservation in the decay, Pyi* = (— 1 )'^^!+^. Angular momentum conservation 
also requires max(J^ — |, 0) < + §• We assume the minimal value of in 

The electromagnetic decay -0 —)■ is described by a term 




(35) 


where AA^ = Ap+ — A^- = ±1, and 0^, 6^ are the azimuthal and polar angles of /i"*" for 
A° (/i“ for A° decays) in the 0 rest frame (see Fig. 16). There are no helicity couplings in 
Eq. ( [ 3 ^ , since they are all equal due to conservation of C and P parities. Therefore, this 
coupling can be set to unity as its magnitude and phase can be absorbed into the other 
helicity couplings which are left free in the ht. The calculation of the 0 decay angles is 


analogous to that of the A* decay angles described above (Eqs. (28)-(30)) 


cos 0^ = - p 
0„ = atan2 


with 


■p». 


(36) 

(-(P,r* X A* 

'i . n M A fhl X {p}\ 

(37) 

A W _ - A*} _ 

"''O — -^0 — 

aP& 

•^3 

(38) 


AS 


and given by Eq. (29 


Collecting terms from the subsequent decays together, the matrix element connecting 
different helicity states of the initial and the final-state particles for the entire A* decay 
chain can be written as 


A/f 

■^'^‘A.o,Ap,AA^, 

b 




X 


E« 


^iXA*4>K ^ 2 


^A* 5 '^•0 


-aJ^.o) (39) 


b 


\A* I -'p 


Terms with different helicities of the initial and hnal-state particles (Ap, AA^) must be 
added incoherently 


M 


/l* |2 


1 + P 




EEh 


(■^^0—+1/2), Ap, AAp 
b 


1 -p^° 


EEh 


(A^O—~l/2), Ap, AAp 
b 


(40) 

where P^° is the A^ polarization, defined as the difference of probabilities for Ayio = -1-1/2 
and —1/2 Em. Eor our choice of the quantization axis for A^ spin, no polarization is 
expected (P^ = 0) from parity conservation in strong interactions which dominate A° 
production at LHCb. 
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4.3 Matrix element for the decay chain 

We now turn to the discussion of —)■ PcjK~, Pcj —)■ ijjp decays, in which we allow more 
than one pentaquark state, j = 1, 2 ,... . Superscripts containing the Pc decay chain name 
without curly brackets, e.g. 0^'=, will denote quantities belonging to this decay chain and 
should not be confused with the superscript “{Pc}” denoting the P^ rest frame, e.g. 0 _ 

With only a few exceptions, we omit the A* decay chain label. 

The weak decay 4° —)■ PcjK~ is described by the term. 




(41) 


where are resonance (he. j) dependent helicity couplings. The helicity of the 

pentaquark state, Ap^, can take values of independently of its spin, Jp^. = |, |,... . 
Therefore, there are two independent helicity couplings to be determined for each Pcj state. 
The above mentioned 0p^, 9P symbols refer to the azimuthal and polar angles of Pc in 

the A^ rest frame (see Fig. 17). 

Similar to Eq. (25), the helicity angle in the Pc decay chain can be calculated as. 


cosOP = 


{lab} .{UO} 


^ = PAi ■ Pp. 


(42) 


The (pp^ angle cannot be set to zero, since we have already dehned the axis in the 
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1/1^ 

rest frame by the (pA* = 0 convention. According to Eq. (18) (xq^ ^ ) we have: 






■ V^bt _ .f, {la.b}x « {lab} 

A*±zo ~ Pa* \Pa* Pao ) Pao ’ 




Xa = 


^A*±zo 

{^?} 
/I* ±20 


a 


The (j)p^ angle can be determined in the rest frame from 


= atan2 x 'PpA) 


The strong decay Pcj —)■ 'tpp is described by a term 


(43) 


(44) 


*1-1 ^<^3 


\p^,\p-xp 


i4>3p,0p^,0)* Rp (M^p), 


(45) 


where 6p^ are the azimuthal and polar angles of the ip in the Pc rest frame (see Fig. 17). 


.... (p1 

They are dehned analogously to Eqs. (28) —(30). The Zq direction is dehned by the 


boost direction from the rest frame, which coincides with the direction. This 

leads to 

(46) 


„„„ fj _ {^c} ^ {±c} 

cos Up^ — —Pp- ■ 


■Pi> 


The azimuthal angle of the ip can now be determined in the Pc rest frame (see Fig. 17) 
from 

= atan2 {-(vP x AA ' p/'l ' PP) ■ (47) 

In Eq. (45), the direction is dehned by the convention that we used in the A^ 


rest frame. Thus, similar to Eq. (29) we have 




? V^bt _ * {lab} / « *v^bl\* vv 

ho±P, - Pao yPAO PPc ) PPc 


{lab} ^ {xlg}N ^ {xl?} 


{Pc} 


a 


XA = 


2:0 ± Pc 


ho±Pc 


(48) 


We have labeled the ip and p helicities, and Ap'=, with the Pc superscript to make it 
clear that the spin quantization axes are different than in the A* decay chain. Since the 
ip is an intermediate particle, this has no consequences after we sum (coherently) over 
A^"" = —1, 0, +1. The proton, however, is a hnal-state particle. Before the Pc terms in the 
matrix element can be added coherently to the A* terms, the A^'= states must be rotated to 
Ap states (dehned in the A* decay chain). The proton helicity axes are different, since the 
proton comes from a decay of different particles in the two decay sequences, the At* and 
Pc- The quantization axes are along the proton direction in the Al* and the Pc rest frames, 
thus antiparallel to the particles recoiling against the proton: the K~ and ip, respectively. 


These directions are preserved when boosting to the proton rest frame (see Fig. 18). Thus, 


32 













Figure 18: Definition of the 9p angle. 


the polar angle between the two proton qnantization axes {6p) can be determined from 
the opening angle between the K~ and ip mesons in the p rest frame, 


Q ^ {p} ^ {P} 

COS Op = ■ p^\ 


(49) 


(A similar problem is discnssed in Ref. 47 


in ^ AT+tt-Yh decays, in the interference of 
RO ^ Z-K+, Z- ^ x,i7r 


where the two different Xci helicity frames 
-P- K*Xci, K* —)■ K^n~ and of 
) The dot product above must be 


contributions, are realigned, 
calculated by operating on the p^ and p^ vectors in the proton rest frame obtained by 
the same sequence of boost transformations, either according to the A* or R, decay chains, 
or even by a direct boost transformation from the lab framej^ 

No azimuthal rotation is needed to align the two proton helicity frames, since the decay 
planes of the A* and the Pc are the same (see Fig. 18). Therefore, the relation between \p 
and Aiy states is 


IV = E A (0.«p.0)'lT'> = 


Op 




(50) 


Thus, the term given by Eq. ( |45| ) must be preceded by 

{Op). 


\Pc -L 1 

Ap —± 2 




AjyAo 


(61) 


Parity conservation in Pcj —)■ ipp decays leads to the following relation 






\^c \r 
A,/. , Ap 


= P^PpPp^^ y^Pj.- 


^Pc 


= Pp 


- 1 ) 




ipp 
Pc ) 


(52) 


^Numerical values of momentum vector components, {px,Py,Pz), depend on the boost sequence taken 
and are related between different boosts via the rotation matrix. However, the dot product between the 
two vectors remains independent of the boost sequences. 
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where Pp^. is the parity of the Pcj state. This relation reduces the number of inde¬ 
pendent helicity couplings to be determined from the data to 2 for Jp^ = ^ and 3 for 
Jpcj ^ Since the helicity couplings enter the matrix element formula as a product, 

the relative magnitude and phase of these two sets must be fixed by a 

Pc , Ap 

convention. For example, can be set to (1,0) for every P^j resonance, in which 


case , i 

APc=+o 


^ develops a meaning of the complex ratio of 

APc — 


2 '^Pc 2 


while all 


^fpc^p^ couplings should have both real and imaginary parts free in the fit. 

A^, >Ap 

The term Rp^.{M^p) describes the ipp invariant mass distribution of the Pcj resonance 
and is given by Eq. (31) after appropriate substitutions. Angular momentum conservation 
limits LJi' in —)• PcjK~ decays to Jp^. ± 1. The angular momentum conservation also 
imposes max(Jp^^, — |, 0) < Lp^^ < Jp^. P |, which is further restricted by the parity 
conservation in the Pcj decays, Pp^ = (—1) . We assume the minimal values of LP 

and of Lp^. in Rp^.{m^p). 

The electromagnetic decay ip —)■ pApT in the Pc decay chain contributes a term 




(53) 


which is the same as Eq. (35), except that since the ip meson comes from the decay of 
different particles in the two decay chains, the azimuthal and polar angle of the muon in 
the Ip rest frame, cpp‘^, 0^'", are different from (pp, 0^ introduced in the A* decay chain. The 
Ip helicity axis is along the boost direction from the Pc to the ip rest frames, which is given 
by 


and so 


cos 


M Pc _ m 
0 — Pv 


qPc _ M . A M 

— Pp Pp. ■ 


The X axis is inherited from the Pc rest frame (Eq. (19)), 


(54) 

( 66 ) 


a 


X 


{A Pc 


20 ap 

{Pc} 


= -p/‘*+(pT“'-pr')p* 


{Pc} c{Pc}.c{Pc} 


= ay - 


PPc} 

PPc} 


(56) 


which leads to 


= atan2 x x^}^^ a, 


M Pc . f. {A 
0 Pfj, 


(57) 


Since the muons are hnal-state particles, their helicity states in the Pc decay chain, 
|A^'=), need to be rotated to the muon helicity states in the A* decay chain, |A^), before 
the Pc matrix element terms can be coherently added to the A* matrix element terms. 
The situation is simpler than for the rotation of the proton helicities discussed above, as 
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le '0 decay in both decay chains. This makes the polar angle 

I 1 

51)) eqnal to zero, which leads to d T (0) = 6.Pc ^ , where 5*, 


the mnons come from t 
(analogous to 9p in Eq. 
is the Kronecker symbol. However, the muon helicity states are not identical since the x 
axes are offset by the azimuthal angle a^. Since the boost to the /i rest frame is the same 
for both decay chains {i.e. always from the 'ijj rest frame), we can determine in the 'ijj 
rest frame 


where and from Eq. (19) 


= atan 2 




X X 


M Pc 


X 


M A* 'm p< 


X 


M A" 


X 


M Pc _ 


= — a 


■ m 


= ~Pp"' + (P: 


M Pc 

Zo±fJ. 

#} 


{hi 


as well as 


X 


M A* _ 


= — a 


M A* 

Zo±fJ. 


^ #} A* 

^Z0±fJ, 


c {?/)} 

= -Pa* 


+ (Pit 


m 


Ffi ) Ffi 5 


p«)pM. 


(58) 

(59) 

(60) 

(61) 

(62) 


The term aligning the muon helicity states between the two reference frames is given by 


ETt.vK.o.o)' = E 


c 0 \ Pc \ — o 


(63) 


The transformation of fi~ states will be similar to that of the states, except that since 
will have the opposite direction, = —a^-. The transformation of |A^+)|A^l) to 
|A^+)|A^-) states will require multiplying the terms for the Pc decay chain by 

(64) 


X^_ol^^ iX^OLjx _ A^)q;^ _ 


i AXn 


An alternative derivation of Eq. (64) is discussed in Ref. (Eqs. (20) —(22) therein) for 
the interference of —)■ K*ip, K* —)■ Ktt and of ZK~, Z —)■ xJjtt {'ip —)■ PP) terms, 
which are analogous to the two decay chains discussed here with the substitution 5° —)■ 

K* A*, Z ^ Pc and vr —)• p. The rotation by about the £~^ direction in the pj rest 


frame in the Z decay chain is incorporated by setting 7 = a^, instead of 7 = 0 in Eq. (53) 
This leads to the same formulae since 


D ^ 




= -DA,aa,(W. 


(65) 


We use the more generic derivation here to demonstrate that the methods of transforming 
the muon and proton helicity states between the two decay chains are the same. 
Collecting terms from the three subsequent decays in the Pc chain together. 


\ r. A 


X 


j 

E« 


A .p p p ^ p p 

V.ApT Ai) A, A ^PcA-A 


Pcn^i^P 7 ’^Pc 
Pc 

ijj ’ '^Pc^'A 


(nPc\nj^<^3-^^P 
AO ) ^ ' p- ' P'' ^ 


(»ft). (66) 
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and adding them coherently to the A* matrix element, via appropriate relation of 
|Ap)|A^+)|A^-) to |Ap'=)|A^+)|A^i:) states as discnssed above, leads to the hnal matrix 
element sqnared 


E E E 




M 


/\0 i 

b 






\^c \ 

Ap , Aj 


AA 

‘^j\0 } Ap , AAp 

b 


(67) 

where we set = 0. As a cross-check, htting the polarization to the data with 
the defanlt A* and model yields a valne consistent with zero, = (—2.0 ± 2.3)% 
(statistical error only). 

Assuming approximate CP symmetry, the helicity couplings for A^ and A° can be made 
equal, but the calculation of the angles requires some care, since parity (P) conservation 
does not change polar (i.e. helicity) angles, but does change azimuthal angles. Thus, not 
only must 'p^+ be used instead of for A° candidates (with and p in the hnal-state) 
in Eqs. (36), (37), (55), (57) and (58), but also all azimuthal angles must be reflected before 

4>k, 4>ii 


entering the matrix element formula: (px 
—)■ —0^'" and a,,, —)■ —a ^ 


-)■ 




-)■ 




-)■ 






It is clear from Eq. (67) that various Ap and Pc resonances interfere in the differential 
distributions. By integrating the matrix element squared over the entire phase space the 
interferences cancel in the integrated rates unless the resonances belong to the same decay 
chain and have the same quantum numbers]^ 


4.4 Reduction of the number of helicity couplings 


A possible reduction of the helicity couplings can be achieved by relating them to the LS 
couplings {Bl,s) using Clebsch-Gordan coefficients 


^As,Ac ~ . 


2L+1 
2 J ^+1 


B 


L,S 


Jb 

As 


Jc 

—Ac 


As — Ac 


X 


L 

0 


As — Ac 


Ja 

As — Ac 


( 68 ) 

and then restricting the L values. Here L is the orbital angular momentum in the decay, 
and S is the total spin of the daughters, S = Jb + Jc (l-^s — Jc\ C S < Jb + Jc)- h fhe 
energy release in the decay, Q = Ma — Mb — Me, is small, Q/Ma <C 1, then higher values 
of L should be suppressed; this effect is usually called “the angular momentum barrier.” 

A* 

Applying this approach to A^ —>■ ipAp decays, the lowest LP value (Amin) corresponds to a 

single possible value of S', thus reducing the number of couplings to fit, from 4 (J/i* = |) 
or 6 (J/i* > |), to just one Bl s coupling per resonance. Accepting also Lcain + 1 values, 
gives three Bb^s couplings to fit per resonance. 


In A° — Pr^K decays, S = Jp^. and = Jp^. ± Taking only the lower 

0 ^b 


'■b ' cj 

value reduces the number of couplings from 2 to 1. 


Since its magnitude and phase 


®For A*i~Pc, the A^io = -1-1/2 interference terms have the opposite effect to the Ay^o = —1/2 interference 
terms. 
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convention can be absorbed into T-i p 

A/, Ap 


bp 

Pc 


'see the discussion in Sec. 


set ij d 

JPcj-2’JPcj 


4.3), one can simply 


= (1, 0) in this approach. 


The reduction of couplings to £t for Pc? —t tjjp decays depends on the spin and parity 

Values of Lp^. must be odd (even) for even 


of the Pcj state. 


S can take values of \ and |. 

1 + 

2 


(odd) Ppcj- For a Jp^. = state, only Lp^^ = 1 is allowed with the two possible values of 

S. Therefore, no reduction of couplings is possible. For a Jp^^ = | state, Lp^. =0,2 are 
allowed, each corresponding to one S value. Therefore, the number of couplings to ht can 
be reduced from 2 to 1 when taking Lp^^ = 0. Gains can be larger for Jp^. > | states. 
Even if no reduction in parameters is achieved, expressing the helicity couplings via 


corresponding Bp^s couplings using Eq. (68) is useful, since it automatically implements the 

parity constraints (Eq. (52)) by restricting possible L values. Since the overall magnitude 

of the matrix element does not affect the normalized signal PDF, and because its overall 

phase also drops out when taking its modulus, we fix the magnitude and phase convention 
yl(1520)J/b 


ylO- 

by setting B \ 

U, 2 


= ( 1 , 0 ). 


5 Details of fitting techniques 


The matrix element given by Eq. ( [67| ), is a 6-dimensional (6D) function of rriKp and 
i7 and depends on the ht parameters, at, which represent independent helicity or LS 
couplings, and possibly masses and widths of resonances, Ai{mKp, f2| at). We perform 
an unbinned maximum likelihood ht of these parameters to the 6D data by minimizing 
—21n£(at) = —2 In 'P(mii'p *,17*1 at) with respect to at’. The signal PDF is obtained 
by multiplying the matrix element squared with the selection efficiency, e{mKp, 17), 




’sig(^”A>. ^1^) ^ \M{mKpyP\'Zt)f ^{mKp)e{mKp,0), 


(69) 


where <F(mxp) is the phase space function equal to pq, where p is the momentum of the 
Kp system (he. A*) in the rest frame, and q is the momentum of K~ in the A* rest 
frame, and /(at) is the normalization integral. 

We use two ht algorithms that were independently coded and that differ in the approach 
used for background subtraction. The sFit procedure uses the sPlot technique described in 
Ref. 31 to subtract background from the log-likelihood sum. It has been used in previous 
LHCb analyses, e.g. measurement of (f)s in B^s decays 48 . The data in the entire 
5480 — 5760 MeV range are passed to the htter. Candidates are assigned “sWeights” (W) 
according to their mj/^px value with the signal and background probabilities determined 
by the ht to the mj/^px distribution. 


— 2 In/ICot) = —2sw hFj \nV{mxpi, l7j| at). 


(70) 


where sw = ^i/ is a constant factor accounting for the effect of the background 

subtraction on the statistical uncertainty. Candidates in the sideband region have negative 
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sWeights which on average compensate for the backgronnd events present in the peak 
region, where events have positive sWeights. From signal simnlations, we see signihcant 
variations of —)■ J/'iljpK~ invariant mass resolntion as fnnctions of cos^^io and cos6j/^. 
The background distributions also vary in intensity and shape with these two variables. No 
strong variation is seen for the other fitting observables. To determine the sWeights, the 
events are divided into 4 | cos9j/^\ x 8 cos6'yio bins. A separate fit to the —)■ J/'ippK~ 
invariant mass distribution of each bin is performed. 

In the sFit approach the total PDF is equal to the signal PDF, as the background is 
subtracted from the log-likelihood sum using sWeights, 

-21n£(at) = - 2sw W lnT’sig(^jrp i, 

i 

= — 2sw E Wi In \M{mKp i, at)|^ -h 2sw In J(at) ^ Wi 

i i 

-2sw y^^Wi\n[^{mKp i)e{mKp i, (71) 

i 

The last term does not depend on the fitted parameters at and is therefore dropped. The 
efficiency still appears in the normalization integral. The integration is done without the 
need to parameterize the efficiency, by summing the matrix element squared over the 
simulated events that are generated uniformly in phase space and passed through the 
detector modeling and the data selection procedure, 

/(at) = J |Al(mi^p, i7| at)|^ <l>(mxp)e(mifp, i7) 

j,f7j|7t')|^, (72) 

j 

where are the weights given to the simulated events to improve the agreement between 
data and simulations. They include particle identification weights obtained from calibration 
samples of A —)■ p7r~ and —)■ K~tt^ as functions of momentum and pseudorapidity of 
the protons and kaons. They also include a weight correcting the overall efficiency for the 
dependence on the charged track multiplicity of events, determined from the ratio of the 
background-subtracted data and the signal simulations. Imperfect description of the A° 
production kinematics in our simulation is corrected in a similar way via a weight that 
depends on the p and px of the A° baryon. The final weights are the ratio of the data and 
the simulations as a function of proton and kaon momenta. 

In the cFit approach, candidates are not weighted (Wi = 1). The data that are fitted 
are restricted to be within a ±2(T mass window around the A^ mass peak, in the interval 
5605.7 < mj/^pK < 5635.8 MeV. The background is represented explicitly in the htted 
PDF, with the integrated background probability set to /3 = 5.4% as determined from the 
fit to the J/il)K~p mass distribution, 

/^(mi^p, i7| at) = (1 -/3)'Psig(mii-p, i7|at) +/?'Pbkg(mxp, f?)- (73) 
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The 6-dimensional backgronnd parameterization V^kgimKp, 1^) is developed nsing the 
background sample in which the candidate invariant mass is more than 5a away from 
the peak, specifically within the intervals 5480.0 — 5583.2 MeV and 5658.3 — 5760.0 MeV. 
We minimize the negative log-likelihood defined as, 


2 In £(75’) = 

-2^1n 

i 

= - 25 ^ 1 n 


(1 _ ^{rriKp i)e{mKp i, ^ ^ '^hkgij^Kp i, ^i) 


\M{mKp i, + 


7(^ 


'bkg 


(1 /^)7bkg i)(^(j^Kp j) 


2iVln J(at) -|- constant. 


(74) 


where N is the number of candidates, and 'P^kgi'^Kp, is the unnormalized background 
density proportional to the density of sideband candidates, with its normalization deter¬ 
mined bjj^ 


/bkg = 


V^kgij^Kp) drriKpdQ oc ^ 


^,MC '^hkgi^Kp j, ^j) 

^ ^{rriKp j)e{mKp 


(75) 


The background term is then efficiency-corrected so it can be added to the efficiency- 
independent signal probability expressed by |W1|^. This way the efficiency parametrization, 
e{mKpi f2), influences only the background component which affects only a tiny part of 
the total PDF, while the efficiency corrections to the signal part enter Eq. (72). 

The efficiency in the background term is assumed to factorize as 


e{mKp, 12) = ei(mKp,cos6');)e2(cos6'^o|mi^p)e3(cos6'j/^|mi^p)e4(0ic|m;^p)e5(0/,|m;^p). 

(76) 

The ei{mxp, cos 9*^ term is obtained from a 2D histogram of the simulated events weighted 
by l/(pg). A bi-cubic interpolation is used to interpolate between bin centers. The other 
terms are again built from 2D histograms, but with each bin divided by the number of 
simulated events in the corresponding rriKp slice to remove the leading dependence on this 
mass, which is already taken into account in the hrst term. The leading variation of the 
efficiency is in the ei(mj^p, cos 6*^) term which is visualized in the normal Dalitz variables 

a) instead of the “rectangular Dalitz plane” variables {rriKp, cos 9\) 




in Fig. 
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used to parameterizeThis variation. 

The background PDF, V^^^^irriKpi 12)/d>(mfCp), is built using the same approach. 


Vkkgi.'^Kp, 12) 
^{niKp) 


Tbkgi (mxp, cos 6';^)Pbkg2(cos 0^0 \mKp) 


X Pbkg3 (cos 0 |//^Xp)Pbkg4 (^iC |//^jrp)Pbkg5 (0p |//^jrp) • 

A visualization of Pbkgi(/^irp 5 cos9\) on the Dalitz plane is shown in Fig. [I^b). 

®Note that the distribution of MC events includes both the ^{rriKp) and e{mKp, P) factors, 
cancel their product in the numerator. 


(77) 

which 
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Figure 19: Parameterized dependence of (a) the relative signal efficiency and of (b) the background 
density on the Dalitz plane. The units of the relative efficiency and of the relative background 
density are arbitrary. 
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